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Abstract. We consider the possibility of observing CP violation in quasi- inclusive 
decays of the type B~ -> K~X, B~ -► K*~X, B° -> K~X and B° -> K*~X, where 
X does not contain strange quarks. We present estimates of rates and asymmetries 
for these decays in the Standard Model and comment on the experimental feasibility 
of observing CP violation in these decays at future B factories. We find the rate 
asymmetries can be quite sizeable. 

I INTRODUCTION 

The possibility of observing large CP violating asymmetries in the decay of B 
mesons motivates the construction of high luminosity B factories at several of the 
world's high energy physics laboratories. The theoretical and the experimental 
signatures of these asymmetries have been extensively discussed elsewhere [1], [2], 
[3], [4], [5]. At asymmetric B factories, it is possible to measure the time dependence 
of B decays and therefore time dependent rate asymmetries of neutral B decays due 
to B — B mixing. The measurement of time dependent asymmetries in the exclusive 
modes B° — > ipK s and B° — > 7r + 7r~ will allow the determination of the angles in 
the Cabbibo-Kobayashi-Maskawa unitarity triangle. This type of CP violation has 
been studied extensively in the literature. 

Another type of CP violation also exists in B decays, direct CP violation in the 
B decay amplitudes. This type of CP violation in B decays has also been discussed 
by several authors although not as extensively. For charged B decays calculation 
of the magnitudes of the effects for some exclusive modes and inclusive modes have 
been carried out [6], [7], [8], [9], [10], [11], [12]. In contrast to asymmetries induced 
by B — B mixing, the magnitudes have large hadronic uncertainties, especially for 
the exclusive modes. Observation of these asymmetries can be used to rule out the 
superweak class of models [13]. 

In this paper we describe several quasi-inclusive experimental signatures which 
could provide useful information on direct CP violation at the high luminosity fa- 
cilities of the future. One of the goals is to increase the number of events available 



at experiments for observing a CP asymmetry. In particular we examine the inclu- 
sive decay of the neutral and the charged B to either a charged K or a charged K* 
meson. By applying the appropriate cut on the kaon (or K*) energy one can isolate 
a signal with little background from b — > c transitions. Furthermore, these quasi- 
inclusive modes are expected to have less hadronic uncertainty than the exclusive 
modes, would have larger branching ratios and, compared to the purely inclusive 
modes they may have larger CP asymmetries. In this paper we will consider modes 
of the type B — > K(K*)X that have the strange quark only in the K(K*) -meson. 

In the sections which follow, we describe the experimental signature and method. 
We then calculate the rates and asymmetries for inclusive B~ — > K~(K*~) and 
B° -> K~(K*-) decays. 

II EXPERIMENTAL SIGNATURES FOR 
QUASI-INCLUSIVE B -> SG* 

In the T(4S') center of mass frame, the momentum of the from quasi-two 

body B decays such as B — > K^~X may have momenta above the kinematic 
limit for mesons from b — > c transitions. This provides an experimental 

signature for b — > sg*, g* — > or g* — > <Zd decays where g* denotes a gluon. This 
kinematic separation between 6 — > c and 6 — > sg* transitions is illustrated by a 
generator level Monte Carlo simulation in Figure 1 for the case of B — > A" * _ . (The 

— > A' - spectrum will be similiar). This experimental signature can be applied 
to the asymmetric energy B factories if one boosts backwards along the z axis into 
the T(AS) center of mass frame. 

Since there is a large background ( "continuum" ) from the non-resonant processes 
e + e~ — > qq where q = u,d, s, c, experimental cuts on the event shape are also 
imposed. To provide additional continuum suppression, the "B reconstruction" 
technique has been employed. The requirement that the kaon and n other pions 
form a system consistent in beam constrained mass and energy with a B meson 
dramatically reduces the background. After these requirements are imposed, one 
searches for an excess in the kaon momentum spectrum above the b — > c region. 
Only one combination per event is chosen. No effort is made to unfold the feed- 
across between submodes with different values of n. 

Methods similar to these have been successfully used by the CLEO II experiment 
to isolate a signal in the inclusive single photon energy spectrum and measure the 
branching fraction for inclusive b — > transitions and to set upper limits on b — > scf) 
transitions [14], [15]. It is clear from these studies that the B reconstruction method 
provides adequate continuum background suppression. 

The decay modes that will be used here are listed below: 

1. B~ -> tfW-TT 

2. B° -> ATW-tt + 
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FIGURE 1. Generated Inclusive B — * K*~ momentum spectrum. The component below 2.0 
GeV/c is due to b — > c decays while the component above 2.0 GeV/c arises from quasi-two body 
b — > sg* decay. The normalization of the & — > c component is reduced by a factor of approximately 
100 so that both components are visible. 



3. B- -> iifW-7r-7r + 

4. 5° -> iv"W-7r+7r 

5. 5° -> iv"W-7r + 7r-7r + 

6. B- -> fSfW-TT + Tr-Tr 

7. B- -> fs:W-7r + 7r-7r + 7r- 

8. 5° -> iv-W-TT + TT-TT + Tr 

In case of multiple entries for a decay mode, we choose the best entry on the 
basis of a x 2 formed from the beam constrained mass and energy difference (i.e. 
X 2 = (Mb/ 5 Mb) 2 + (AE /5AE) 2 ). In case of multiple decay modes per event, the 
best decay mode candidate is picked on the basis of the same x 2 - 

Cross-feed between different b — > sg decay modes (i.e. the misclassification of 
decay modes) provided the is correctly identified, is not a concern as the goal 

is to extract an inclusive signal. The purpose of the B reconstruction method is 
to reduce continuum background. As the multiplicity of the decay mode increases, 
however, the probability of misrecontruction will increase. 

The signal is isolated as excess production in the high momentum signal 

region (2.0 < p K (*) < 2.7 GeV) above continuum background. To reduce con- 
tamination from high momentum B — > n~(p~) production and residual b — > c 



background, we assume the presence of a high momentum particle identification 
system as will be employed in the BABAR, BELLE, and CLEO III experiments. 



We propose to measure the asymmetry N(K^ + —K^~)/N(K ( -*" I+ +K ( -*^) where 
j((*)± originates from a partially reconstructed B decay such as B — > K^~(nn) 
where the additional pions have net charge and n < 4 and one neutral pion is 
allowed and 2.7 > p(K^~) > 2.0 GeV. We assume that the contribution from 
B — > K~r\ X decays has been removed by cutting on the i] region in X mass. It is 
possible that the anomalously large rate from this source [16,17] could dilute the 
asymmetry. 
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FIGURE 2. Monte Carlo simulation of inclusive B — > K*~X signal with the B reconstruction 
method: (a) The beam constrained mass distribution (b) The distribution of energy difference 
(c) The K~ir° invariant mass after selecting on energy difference and beam constrained mass 
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III EFFECTIVE HAMILTONIAN 

In the Standard Model (SM) the amplitudes for hadronic B decays of the type 
b — > qff are generated by the following effective Hamiltonian [18]: 

Q 10 

H q eff = ^[v fb v; q ( Cl oi f + C2 ot f ) - J2(v ub v: qC ? + v cb v; q c< + v tb v; q 4)oi} + h.c. , (i) 

where the superscript u, c, t indicates the internal quark, / can be u or c quark. 
q can be either a d or a s quark depending on whether the decay is a AS* = or 
AS = — 1 process. The operators 0\ are defined as 

0% = q al „Lf p f^Lb a , 0\ f = q^Lff-fLb , 

3 % = q^Lbq'^L(R)q' , 4 % = q^.Lb^^R)^ , (2) 

3 3 
Cr,9 = 7pl»Lbe ql q'^R{L)q' , Og jl0 = -q^^Lb^q'^^RiL)^ , 

where R(L) = 1 ± 75, and q' is summed over u, d, and s. 0\ are the tree level and 
QCD corrected operators. are the strong gluon induced penguin operators, 

and operators O 7 _i are due to 7 and Z exchange (electroweak penguins), and 
"box" diagrams at loop level. The Wilson coefficients c{ are defined at the scale 
pL m b and have been evaluated to next-to-leading order in QCD. The c\ are the 
regularization scheme independent values obtained in Ref. [9]. We give the non-zero 
c{ below for m t = 176 GeV, a s (mz) = 0.117, and /x = m b = 5 GeV, 

Cl = -0.307 , c 2 = 1.147 , 4 = 0.017 , c\ = -0.037 , c* = 0.010 , c\ = -0.045 , 
4 = -1.24 x 10~ 5 , 4 = 3 - 77 x 10 ~ 4 , 4 = -0.010 , 4o = 2 - 06 x 10 " 3 , 

4; 5 C = -c u 4/n c = p:> c /n c , <%$ = pt, cZ w = o (3) 



where N c is the number of color. The leading contributions to P l se are given by: 

^)(iV cCl + c 2 )(f 



PI = (^hff+GKft? 2 )) and Pt = (^)(iV cCl + c 2 )(f + G( mi ,^q 2 )). The 



function G(m, //, g 2 ) is given by 



G(m, n,q 2 ) = 4 x(l - a) In — dx . (4) 

•/ 



All the above coefficients are obtained up to one loop order in electroweak inter- 
actions. The momentum g is the momentum carried by the virtual gluon in the 
penguin diagram. When q 2 > 4m 2 , G(m,fj,,q 2 ) becomes imaginary. In our cal- 
culation, we use m u = 5 MeV, rrid = 7 MeV, m s = 200 MeV, m c = 1.35 GeV 
[19,20]. 

We assume that the final state phases calculated at the quark level will be a good 
approximation to the sizes and the signs of the FSI phases at the hadronic level 
for quasi-inclusive decays when the final state particles are quite energetic as is the 
case for the B decays in the kinematic range of experimental interest [6]. 



IV MATRIX ELEMENTS FOR B -> K X AND 

We proceed to calculate the matrix elements of the form < KX\H e ff\B > which 
represents the process B — > ifX and where i? e // has been described above. The 
effective Hamiltonian consists of operators with a current x current structure. 
Pairs of such operators can be expressed in terms of color singlet and color octet 
structures which lead to color singlet and color octet matrix elements. In the 
factorization approximation, one separates out the currents in the operators by 
inserting the vacuum state and neglecting any QCD interactions between the two 
currents. The basis for this approximation is that, if the quark pair created by one of 
the currents carries large energy then it will not have significant QCD interactions. 
In this approximation the color octet matrix element does not contribute because 
it cannot be expressed in a factorizable color singlet form. In our case, since the 
energy of the quark pairs that either creates the K or the X state is rather large, 
factorization is likely to be a good first approximation. To accommodate some 
deviation from this approximation we treat N c , the number of colors that enter in 
the calculation of the matrix elements, as a free parameter. In our calculation we 
will see how our results vary with different choices of iV c . The value of N c ~ 2 
is suggested by experimental data on low multiplicity hadronic B decays [3]. The 
amplitude for B — > K^X can in general be split into a three body and a two 
body part. Detailed expressions for the matrix elements, decay distributions and 
asymmetries can be found in [21] 

V RESULTS AND DISCUSSION 

In this section we discuss the results of our calculations. We find that there 
can be significant asymmetries in B — > K(K*)X decays especially in the region 
Ek > 2 GeV which is also the region where an experimental signal for such decays 
can be isolated. The branching ratios are of order O(10" 4 ) which are within reach 
for future B factories. The contribution of the amplitude with the top quark in 
the loop accounts for 60-75% of the inclusive branching fraction. However, since 
the top quark amplitude is large and has no absorptive part in contrast to the c 
quark amplitude, the top quark contribution reduces the net CP asymmetry from 
30-50% to about 10%. This calculation includes the contribution from electroweak 
penguins. We find that the electroweak penguin contributions increase the decay 
rates by 10-20% but reduce the overall asymmetry by 20-30%. The main sources 
of uncertainties in our calculation are discussed extensively in [21]. 

The asymmetries are sensitive to the values of the Wolfenstein parameters p and 
T). The existing constraints on the values of p and r\ come from measurements of 
|Kt6|/|K*|) e K in the K system and AM^. (See Ref. [22] for a recent review). In 
our calculation we will use = 170 MeV and choose (p = —0.15, 77 = 0.33). 
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FIGURE 3. Predicted Asymmetries for B~ —* K ~X and B~ — > K*~X as a function of the 
kaon energy. The three sets of curves indicate the sensitivity of the asymmetry to the value of 
N c . The values N c — 2, 3, oo are considered. 



TABLE 1. Integrated decay rates and asymmetries for B — > K^X Decay 



Process Branching Ratio (1.65 x 10 4 ) Integrated Asymmetry 

B~^K~X 1.02,0.79,1.20 -0.10,-0.11,-0.050 

B~ -» K-X{E K > 2.1GcV) 0.81,0.74,0.77 -0.12,-0.12,-0.07 

Tf^K-X 0.6,0.7,0.8 -0.12,-0.12,-0.13 

B-^K*-X 1.37,1.24,2.30 -0.11,-0.14,-0.11 

B~ -» K*-X(E K , > 2.1GeV) 1.05,1.16,1.67 -0.14,-0.15,-0.14 

B°->K*-X 1.05,1.16,1.39 -0.15,-0.15,-0.16 



In Fig. 3 we show the asymmetries for K and K* in the final state in charged B 
decays for different values of N c . Variation of the asymmetries with the different 
inputs in our calculation are presented in detail in [21]. 

In Table. 1 we give the branching fractions and the integrated asymmetries for 
the inclusive decays for different N c , q 2 = ml/2 ( q is the gluon momentum in the 
two body part of the amplitude ), fs = 170 MeV, p = —0.15, 77 = 0.33. For the 
charged B decays we also show the decay rates and asymmetries for Ek > 2 (2.1) 
GeV as that is the region of the signal. 

The above figures show that there can be significant asymmetries in B — > K^X 
decays, especially in the region E K > 2 GeV which is the region of experimental 
sensitivity for such decays. As already mentioned, our calculation is not free of the- 
oretical uncertainties. Two strong assumptions used in our calculation are the use 



of quark level strong phases for the FSI phases at the hadronic level and the choice 
of the value of the gluon momentum q 2 in the two body decays. Other uncertainties 
from the use of different heavy to light form factors, the use of factorization, the 
model of the B meson wavefunction, the value of the charm quark mass and the 
choice of the renormalization scale \x have smaller effects on the asymmetries [21]. 

VI CONCLUSION 

We find significant direct CP violation in the inclusive decay B — > K~X and 
B — > K*~X for 2.7 > E K (*) > 2.0 GeV. The branching fractions are in the 10~ 4 
range and the CP asymmetries may be sizeable. These asymmetries should be 
observable at future B factories and could be used to rule out the superweak class 
of models. 
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